Analytical theory of the dressed bound state in highly polarized Fermi gases 
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We present an analytical treatment of a single J. atom within a Fermi sea of j atoms, when the 
interaction is strong enough to produce a bound state, dressed by the Fermi sea. Our method makes 
use of a diagrammatic analysis, with the involved diagrams taking only into account at most two 
particle-hole pairs excitations. The agreement with existing Monte-Carlo results is excellent. In 
the BEC limit our equation reduces exactly to the Skorniakov and Ter-Martirosian equation. We 
present results when f and J. atoms have different masses, which is of interest for experiments in 
progress. 
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The field of ultracold Fermi gases [H has turned re- 
cently toward very interesting and unexplored physical 
domains for superfluidity. Indeed it is possible to have 
stable mixtures of two fermionic species with different 
particle numbers Q (called polarized gases) and also dif- 
ferent masses. When quantum degeneracy is reached, 
these fermions may form Cooper pairs, leading to BCS- 
like superfluidity. In the regime where these pairs are 
very tightly bound, forming essentially small molecules, 
one obtains the Bose-Einstein condensation (BEC) of 
these resulting bosons. However when the attractive in- 
teraction gets weaker, which can be achieved at will ex- 
perimentally through Feshbach resonances [l| , these pairs 
are overlaping increasingly and one has to explore the 
whole extent of the BEC-BCS crossover. Polarization 
has been shown to be detrimental to superfluidity in this 
crossover, but one could also hope to discover new super- 
fluid phases, like the FFLO phases, although experimen- 
tal results are up to now negative in this respect. 

A particularly attractive limiting case is the one of very 
strong polarization, which is at the same time easier to 
handle but provides also the possibility to understand 
quantitatively the physics at lower polarization 0, 0, S] • 
This is the situation where a single fermion (say j-atom) 
with mass mi is in the presence of a Fermi sea of an- 
other, say 1, fermion species with mass m = mt, Fermi 
wave vector kp and scattering length a between ] and 
I atoms. In the absence of bound state, we have shown 
recently @ that the exact solution of this problem is ob- 
tained as the limit of an extremely rapidly convergent se- 
ries of approximations. These successive approximations 
amount to restrict the Hilbert space of the excited states 
of the system to one, two,...,n,.. particle-hole pairs. In 
practice the first order approximation, which coincides 
with the standard ladder approximation is already quite 
satisfactory. The second order one, where at most two 
particle-hole pairs are coming in, gives an essentially ex- 
act answer, as it may be confirmed by comparison with 
QMC calculations, when available. 

In this paper we address analytically, by a diagram- 



matic extension of the above treatment, the case where 
the attractive interaction is strong enough to lead to a 
bound state, which is merely a molecule in the 'BEC 
limit of very strong attraction hp a — > 0+. This regime is 
of essential importance [B| for the understanding of the 
whole phase diagram, since when the density of j atoms 
increases, the corresponding bound states form a Bose- 
Einstein condensate. We will calculate in this regime the 
chemical potential and the effective mass, with results for 
the case of equal masses mj = TO| in remarkable agree- 
ment with known results from QMC calculations @] • 
This allows us to provide accurate answers in the case of 
very high current interest where the mass ratio r = m^/m 
is different from 1. In the BEC limit our results are exact, 
since they reduce to the Skorniakov and Ter-Martirosian 
[§] equation. 

It is physically quite clear [9| that the existence of a 
bound state corresponds to a pole in the vertex corre- 
sponding to the forward scattering of the |-atom and an 
"f-atom. Precisely one can show that such a pole gives 
rise to a singularity in the self-energy £1 (p, ui) for lu < 0, 
leading to a non-zero density ^ [lOj . Since we have 
a single .[-atom iVj = 1, we are, in the thermodynamic 
limit, just at the border between a zero density = 
and a non-zero density ^ 0. Considering the evolu- 
tion of the singularities of £j.(p, w ) when fi± is increasing, 
starting from very large negative values, we find that, for 
N± = La singularity has just reached the border w = 0, 
as in [llj. Moreover when the bound state just appears, 
this corresponds to a zero total energy = for the two 
atoms. Hence the energies of the two scattering atoms are 
zero in our case of interest. The physical fi± is the low- 
est value either giving rise to a bound state or satisfying 
the condition /ij = £(0,0), used in [ill ], when no bound 
state exists. With respect to the momentum variables it 
is clear physically that the lowest energy for the pole is 
obtained when the total momentum q of the two atoms 
is zero. More precisely the dependence of the chemical 
potential on this momentum will give the effective mass, 
which is found positive in the physical range. 
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FIG . 1 : Diagrammatic representation of 7 (k, q) and 
ff(k,q,q') 

Our approximation is the following: in all the diagrams 
we draw, we have at most two explicit propagator lines 
corresponding to "f-atoms running forward through the 
diagram. Naturally we have also the propagator line of 
our single |-atom which runs forward throughout the di- 
agram. We could generalize to 3, 4 • • • explicit prop- 
agators, this series of approximations converging very 
rapidly to the exact result, just as in @, but again the 
present approximation is quite enough. On completely 
general grounds all the diagrams we can draw for the ver- 
tex begin with an elementary interaction between our |- 
atom and an T-atom. However, since we let the strength g 
of this interaction go to zero g — > 0, we have to repeat this 
scattering an infinite number of times, any finite number 
giving a zero contribution. The summation of this series 
gives a factor T 2 (q, il), instead of the factor g we would 
have for a single interaction. It depends only on the total 
momentum q and energy O of the two scattering atoms, 
and not on the variables of each atom separately. In vac- 
uum we would have explicitly T 2 (q, fi) = (q, f2) = 

r -i -1 

(27r/m r ) a -1 — y/2m r (q 2 /2M — Q) (i.e. essentially 
the scattering amplitude of the two atoms with total mass 
M and reduced mass m r ), but here we have to calculate 
this quantity in the presence of the Fermi sea, with chem- 
ical potentials fi-\ and /i^. In a first level approximation, 
where we would have at most a single propagator line for 
|-atoms, Ti would be the only contribution. However at 
our second level approximation, with two f-atoms prop- 
agators, we have also to consider processes where, as a 
first interaction or after T 2 , our single j-atom interacts 
with another "[-atom. 

Before proceeding let us consider first what happens 
when we calculate within our approximation the self- 
energy S| (0, 0) relevant for the case without bound state 
In this case the essential difference is that, in order 
to obtain E from the vertex [Io| , we should not set il = 
but rather close all the diagrams by a j propagator run- 
ning backward. This gives an additional factor Goj (q, 
and we have then to sum over q and Q. We sum over f2 
by closing the f2 contour (which runs on the imaginary 
frequency axis) at infinity in the half-plane Re f2 < 0. 
We do similarly for the other f propagator going back- 
ward. On the other hand, for all the explicit intermediate 
t propagators going forward (that is in the same direc- 
tion as our entering J, propagator), we close the contour 



in the Re a; > half-plane. It can then be seen that, due 
to the explicit f and J. intermediate propagators, all the 
frequency integrations lead merely to on-the-shell eval- 
uations for the remaining factors. This is naturally so 
only when these propagators give rise to poles in the cor- 
responding frequency integration domain. Otherwise the 
result is zero. This leads to the constraint i > fcj? on the 
wave vectors for the forward propagators, and q < kp for 
the ones going backward. In this way we have been able 
to rederive exactly the equation ruling the self-energy in 
6] and accordingly all the results of p for 

Coming back to our present problem, the only differ- 
ence is that, instead of having the on-the-shell evaluation 
il = e q — Hi, with e q = q 2 /2m, we have merely to set 
O = 0. In the following we will have this implicit un- 
derstanding for the variable named q. Otherwise, by 
proceeding in the same way as above, we will have on- 
the-shell evaluations for all the other frequency variables 
entering our diagrams, with the same constraints ki > kp 
and qi < kp on the momenta. Hence, since all the fre- 
quencies are determined in this way, we refer in the fol- 
lowing only to the momentum variables. 

We need to write a general equation for the vertex 
7(k, q) with entering momenta k for the f-atom and q— k 
for the |-atom. Except for a factor, this quantity is ac- 
tually quite similar to ak q in [6], the essential difference 
being the f2 = 0, instead of on-the-shell, value for the q 
variable. A first contribution is naturally T2(q, 0). How- 
ever we have also the possibility that the first interaction 
of the j-atom is with another f-atom, with momentum q'. 
For fixed value of q', we call this contribution H(k, q, q'). 
Since q' is an internal variable, one has to sum over it 
to obtain the contribution to 7(k, q). Finally, after the 
repeated scattering described by T 2 (q, 0), we have the 
possibility to have the processes described by H. These 
two parts will be linked by an f and a j propagator. In- 
tegration over the frequency of the f propagator gives a 
factor — 1/-E£q coming from the on-the-shell evaluation 
of the J, one. In this way we obtain: 



7(k,q) = T a (q,0) 



with El 



1 + 



E 

k'q' 



ff(k',q,q') 



E, 



(i) 

k'q 



-5>(k,q,q')(l) 



Ek-q + £k - Mt and E p = p 2 /2m l , 

(i) 



-kq - I Ml 

quite analogous to E££ in 

We want to find under which conditions 7(k, 0) di- 
verges. This can naturally occur if [T 2 (0, 0)] _1 = 0, but 
this is just the first level approximation. However this 
can also arise from a divergence of H, and we will look 
for the equation corresponding to this condition. Quite 
generally we have for 7 and H a set of coupled linear 
equations. The divergence is obtained when there is a 
solution for the homogeneous part of these equations. In 
the following we retain only the terms contributing to 
this homogeneous part and we omit the other ones. 
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We have now to write for H (k, q, q') an equation analo- 
gous to Eq. U) . However, since we accept only two | prop- 
agators, we do not have to go to higher order vertices, and 
our equation will be closed. First of all, the l-atom and f- 
atom, having a first interaction, will interact repeatedly, 
giving rise to a common factor T 2 (q + q' k, e q / — ek). 
Then, after this process, we may have another H vertex, 
giving rise to a term completely analogous to the second 
one in Eq.©. However we have also the possibility that, 
after this repeated interaction, the involved f-atom re- 
combines with a hole (this corresponds to closing the T 2 
by a backward q' propagator). Afterwards the q' T-atom 
has disappeared and all the possible remaining processes 
are described by 7(k, q). There is again a factor 
coming from the propagator linking these two parts. Fi- 
nally when we replace 7(k, q) by its explicit expression 
from Eq.©, the first term does not contain H and we 
do not retain it. In this way, keeping only the homo- 
geneous part of the equation, we obtain that the bound 
state appears when: 



■e k )\ 1 iJ(k,q,q') = 



[T 2 (q + q'-k, e? , 

~W) Ji}) E q "#(k,q,q ) 



k' ^kk'qq' 



- w ^q v , ~ * 
^kq 

?2(q,0) V ff(k',q,q") 

p(i) 2^ p(i) 

^kq k'q" -^k'q 



(2) 



r>(2) 

J kk'qq' 



with EiZ_, = + Sk+k'-q-q' + £k + £k' - £q' - A 1 ? ■ In 



this equation q appears as a parameter and again the low- 
est energy bound state is obtained when q = 0. Setting 
H(k, 0, q') = Ho(k, q') we find the following equation for 
the appearance of the lowest energy bound state: 



PMq- 



k, e q 



a,,)]- 1 if (k,q) 



k' 



,(2) 

kk'Oq 



E, 



(3) 



1 Vff (h n>\ I ^(0,0) V go(k',qQ 
_^tf (k,q) + — ^ 

kO q' ^kO k'q' C k'0 



E, 



We show now analytically that, in the BEC limit, our 
equation leads systematically to the exact result, what- 
ever the mass ratio r. 

In this BEC limit our equation simplifies since it can 
be seen as the limit kp — ► at fixed a. Since by definition 
q < fcp, we are left with a function Ho(k, 0) = h(k) of the 
single variable k (actually only of its modulus for symme- 
try reasons) . Then the summation over q' in the last two 
terms of Eq.© gives merely a factor n-\ = k F /(6ir 2 ). On 
the other hand, when we compare the first two terms in 
the right-hand side of Eq.©, the only difference is that 
the typical range for the summation over k' is k' ~ 1/a 
(as it can be seen explicitly below) while it is q' ~ kp 
for the summation over q', as we have just seen. Hence 
the second term is of order (kpa) 3 compared to the first, 
and accordingly completely negligible in the BEC limit 



(actually we find numerically that, even for the lowest 
relevant values of X/kpa this term gives a very small con- 
tribution). Accordingly we are left with: 



[T 2 (k-e fe )]- 1 /i(k): 



fe(k') fc| r 2 (O,0) ^fo(k') 

k' Ir'f 



^ r( 2 ) 6tt 2 e {1) 

k' -^kk'OO ^kO 



v k'0 



We are looking for a solution with the form 
P = \i i [\/Ef = £b/Ep + 1 + P(m/m r )kpa which is 
an expansion in powers of kpa, with e?, = l/2m r a 2 the 
molecule binding energy and Ep = k F /2m the Fermi 
energy. Taking into account the general expression 



of T 2 (k,w) which gives [m r T 2 (0, 0)/(2na)]~ 



1 



(2kpa/ir) 1 + y/ (p — l)m r /marctan(Y / (p — l)m r /m) 
and substituting the expansion for p gives the expan- 
sion [r7i r T 2 (O,0)/(27ra)] _1 = -(k F a) 3 [f3/2 + 2/(3tt)]. 
The lowest level approximation corresponds to write 
[T 2 (0, CI)] -1 = 0. This leads to (3 = -4/(3vr) which corre- 
sponds to the Born approximation for the f-atom-dimcr 
scattering length. 

While in the BEC limit we have, in the last term of 
Eq.©, to evaluate T 2 (0,0) carefully since it is large, 
all the other terms can be evaluated to lowest order. 



(m/m r )/ (kpa) 2 in E^ 



(i) 



kO 



Hence we set \pi\/Ep 
\ tH \ + k 2 /2m r andin^k'oo = \pi\ + {k+k') 2 /2m L + (k 2 + 
k' 2 )/2m. Similarly we have [m r T 2 (k, — ek)/(2iTa)] = 
1 — yl + R(ka) 2 , where R = Tn r /rn^ is the ratio of m r 
to the reduced mass = m(m + to j ) /(2m + m i ) of the 
f-atom-dimer system (R = 3/4 for equal masses). We 
take the reduced variable x = ak and make the change 
of function h(k) = C f(x)/x 2 where C is a constant to 
be determined just below. We find from Eq.©: 



irR 1 + VT 



Rx 2 } 1 f(x) = 



(5) 



1 

'2k Jo 



dV 1 dny l + x 2 + y 2 + R>x.y 



1 



4 + 3tt/3 1 + x 2 



dy 



m 

i + y 2 



where R' = 2m r /m^ and the angular average over y 
can be easily performed. We can choose C so that 
/(0) = -3ir/3/(2R). Then, writing Eq.© for x = 0, we 
find that the coefficient of 1/(1 + x 2 ) (including the inte- 
gral) in the last term of Eq.© is just equal to tt. Hence 
Eq. © coincides exactly with the equation, suitably gen- 
eralized to the case of unequal masses, found by Skorni- 
akov and Ter-Martirosian [8| for the scattering amplitude 
to obtain the fcrmion-dimcr scattering length a^. Solving 
Eq.© provides the exact result for 03 = a/(0). From the 
above relation this gives — ~2Ra^/(37ra). When this 
is inserted in the above expression for \pi\, we obtain a 
contribution —2TTn^a^/mJ from the fermion-dimer scat- 
tering, which is precisely the one resulting from a mean- 
field argument, exact in this limit. Hence our equation 
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FIG. 2: Reduced chemical potential a = (/SnmJ /2mEpkFa) 
(tb + Ef — | A* J. I) as a function of 1/hFO, for various mass ratios. 

provides the exact answer in the BEC limit by reducing 
to the Skorniakov and Ter-Martirosian equation. 

We have solved numerically Eq.® in the general case. 
The results for m are given in Fig. [2 or rather we plot 
a = {inmj /2mEFkFa) (e^ + EF — \ni |) which reduces to 
03/0 in the BEC limit and allows to magnify small differ- 
ences. For equal masses mi — m we find that, basically 
down to 1/kpa ~ 2, a is essentially constant, almost 
equal to its BEC value 03/0 = 1.18 (we display also the 
first level result [T 2 (0,0)] _1 = which goes to the Born 
result a^/a — 8/3 in the BEC limit). This is exactly 
what is found by QMC calculations 0, 0|- More pre- 
cisely we find the actual value of a is slightly higher, the 
decrease toward the BEC limit being quite slow, behav- 
ing as kpa. We also find that the bound state appears for 
1/fcpa ~ 0.88 in perfect agreement with QMC (all the 
more remarkable since the two chemical potential curves 
for the polaron and the bound state cross at a very small 
angle, which makes this value quite sensitive). Another 
very striking check of the precision of our results is found 
for mj = 00, where the exact result is known [ill ]. For 
most of the range the agreement is within a few 10~ 3 
and the difference is barely seen even in our blown up 
Fig. O A very interesting feature of this limit is that the 
convergence toward the BEC result a = 1 is as {kpa) 2 , 
faster than in the general case. This allows to understand 
qualitatively why, for the equal mass case, which is not 
much different, one obtains also a fairly slow variation. 
The case of physical interest r = 6.64 for the 40 K- 6 Li 
mixture is also shown and is actually very close to this 
limit. We display also the r = 1/6.64 ~ 0.15 result. 
Note finally that, in agreement with [f|, the quite sim- 
ple q = approximation of Eq. ([3]) gives results in fairly 
good agreement with our exact numerical treatment. 

Finally we show in Fig. [3] our results for the effective 
mass for equal masses = ra. The agreement with di- 
agrammatic QMC results Q is again quite remarkable, 
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FIG. 3: Reduced inverse effective mass rriy/m* of the polaron 
and of the bound state as a function of 1/kpa, at the first 
and second level approximations. The dashed lines indicate 
the location of the appearance of the bound state, both at the 
first level 1/kpa ~ 1.27 and the second level 1/kpa ~ 0.88 
approximations. 

both on the polaron side and on the bound state side. It 
is noteworthy that the effective mass is essentially con- 
tinuous at the transition. It is tempting to speculate that 
this is an exact result, which could be checked experimen- 
tally. This is physically reasonable since, at the threshold 
for bound state appearance, the polaron and the bound 
state are physically identical objects. This physical fact 
is supported by the to^ = 00 explicit case. Note also that 
when the effective mass is calculated beyond this transi- 
tion mi/m* becomes negative, as seen in Fig. [3J which 
signals an instability. As expected physically the transi- 
tion between polaron and bound state occurs before the 
occurrence of this instability. For — 00 the bound 
state appears naturally for l/kpa = 0. On the other 
hand when one gets to lighter m^, the transition value 
for l/kpa goes toward higher and higher values. This is 
already seen easily at our first level approximation which, 
although not accurate as seen in Fig. [3j are nevertheless 
qualitatively correct. 

In conclusion we have shown that our equations pro- 
vide an essentially exact analytical description of the 
molecular state dressed by a Fermi sea. 
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